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Abstract 



We construct the Continuous Wavelet Transform (CWT) on the homogeneous 
space (Cartan domain) D4 = 50(4, 2) /{SO (4) x 50(2)) of the conformal group 
■ 50(4,2) (locally isomorphic to SU{2,2)) in 1+3 dimensions. The manifold B4 

can be mapped one-to-one onto the future tube domain of the complex 
Minkowski space through a Cayley transformation, where other kind of (elec- 
^T) ■ tromagnetic) wavelets have already been proposed in the literature. We study 

, the unitary irreducible representations of the conformal group on the Hilbert 

spaces L 2 t (D4, dv\) and L^{C^_,di>x) of square integrable holomorphic functions 
with scale dimension A and continuous mass spectrum, prove the isomorphism 
(equivariance) between both Hilbert spaces, admissibility and tight-frame con- 
ditions, provide reconstruction formulas and orthonormal basis of homogeneous 
^ ■ polynomials and discuss symmetry properties and the Euclidean limit of the 

proposed conformal wavelets. For that purpose, we firstly state and prove a A- 
extension of Schwinger's Master Theorem (SMT), which turns out to be a useful 
mathematical tool for us, particularly as a generating function for the unitary- 
representation functions of the conformal group and for the derivation of the 
reproducing (Bergman) kernel of L 2 J (B 4 , du x ). SMT is related to MacMahon's 
Master Theorem (MMT) and an extension of both in terms of Louck's SU{N) 
solid harmonics is also provided for completeness. Convergence conditions are 
also studied. 
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1 Introduction 



Since the pioneer work of Grossmann, Morlet and Paul PQ, several extensions of the 
standard Continuous Wavelet Transform (CWT) on M. (traditionally based on the affine 
group of time translations and dilations, see e.g. [2j [3]) to general manifolds X have 
been constructed (see e.g. [U |5] for general reviews and [6j [7] for recent papers on WT 
and Gabor systems on homogeneous manifolds). Particular interesting examples are the 
construction of CWT on spheres S^ -1 , by means of an appropriate unitary representation 
of the Lorentz group in N + 1 dimensions SO{N, 1) El EI], and on the upper sheet 
of the two-sheeted hyperboloid H 2 [11], or its stereographical projection onto the open 
unit disk 

Di = 50(1, 2)/SO{2) = SU(1, 1)/U(1). (1) 

The basic ingredient in all these constructions is a group of transformations G which 
contains dilations and motions on X, together with a transitive action of G on X. 

In this article we shall consider the 15-parameter conformal group G = 50(4,2) 
in 1+3 dimensions and its natural action on the Minkowski spacetime. The fact that 
the conformal group contains space-time dilations and translations leads to a natural 
generalization of the standard CWT for signals on the real line to a higher dimensional 
manifold. Actually, the conformal group 50(4, 2) consists of Poincare transformations 
(space-time translations and Lorentz relativistic rotations and boosts), augmented by 
dilations and relativistic uniform accelerations, which can also be seen as a sort of local 
(point-dependent) scale transformations (see later on Section EJ). 

The conformal group 50(4,2) (or its four-covering SU {2,2)) has been recognized 
as a symmetry of Maxwell theory of electromagnetism without sources since [121 fT3] . 
Electromagnetic waves turn out to be written as superpositions of a particular set of 
conformal wavelets [THESE!]. Thus, conformal wavelets provide a local spacetime-scale 
analysis of electromagnetic waves in much the same way as standard wavelets provide a 
time-scale analysis of time signals. In these works, electromagnetic waves are analytically 
continued or extended from real to complex spacetime and they are obtained from a single 
mother wavelet by applying conformal transformations of space and time. 

Here we shall deal with a different type of conformal wavelets, although we shall work 
in complex spacetime too. Besides the above massless representations of 50(4,2) on the 
electromagnetic field, the conformal group has other representations with continuous mass 
spectrum labelled by the representations of the stability subgroup 50(4) x 50(2): the 
two spins Si,s 2 G N/2 and the scale dimension A G N of the corresponding field [T7] . 
We shall restrict ourselves to scalar fields (si = s 2 = 0) for the sake of simplicity. After 
a reminder of these representations, we provide admissibility conditions, tight wavelet 
frames and reconstruction formulas for functions on the complex Cartan domain or Lie 
ball (see [T8J for a general discussion on these classical complex domains) 

© 4 = 50(4, 2)/(50(4) x 50(2)) = SU{2, 2)/S{U{2) x 17(2)), (2) 

which is the four-dimensional analogue of the open unit disk Di abovementioned. This 
domain can be mapped one-to-one onto the forward/future tube domain (the four- 
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dimensional analogue of the Poincare/Lobachevsky/hyperbolic upper half-plane C+) of 
the complex Minkowski space through a Cayley transformation. For completeness, we 
also provide an isometric (equivariant) map between the Hilbert spaces of holomorphic 
functions on O4 and C+, where we enjoy more physical intuition. 

In order to prove admissibility conditions and reconstruction formulas, an extension 
of the traditional Schwinger's Master Theorem (SMT) [19] will show up as a useful math- 
ematical tool for us. Schwinger's inner product formula turns out to be essentially equiv- 
alent to MacMahon's Master Theorem (MMT) [20J, which is one of the fundamental 
results in combinatorial analysis. A quantum analogue of the MMT has also been con- 
structed [21J and related to a quantum generalization of the boson-fermion correspondence 
of Physics. Moreover, an extension of the classical MMT [20J was proved in [22J by using 
the permutation group. The unification of SMT and MMT into a single form by using 
properties of the so-called SU(N) solid harmonics [23j|2U[25] (a generalization of Wigner's 
D-matrices for 577(2), see e.g. [25]), was pointed out by Louck in |23j. The combined 
MacMahon-Schwinger's Master Theorem provides a generating function for the diagonal 
elements, the trace, and the representation functions of the so-called totally symmetric 
unitary representations of the compact unitary group U(N) j23j |2H 125] . 

In this article we shall state and prove a A-extension of the SMT by using the above- 
mentioned SU(N) solid harmonics of [24| [25] . This A-extension of the SMT will appear 
to be useful as a generating function for the unitary-representation functions of the non- 
compact special pseudo-unitary group 577 (N, N) and for the computation of the repro- 
ducing (Bergman) kernel. We shall concentrate on the N = 2 case, i.e., on the conformal 
group 577(2, 2) (the general case N > 2 is discussed in the Appendix |A|) . which will 
be essential in the development of conformal wavelets for fields with continuum mass 
spectrum. 

The paper is organized as follows. In Section [2] we remind Schwinger and MacMahon's 
Master Theorems and state and prove a A-extension of Schwinger's formula. The general- 
ization to matrices X of size N > 2 is also discussed for completeness in the Appendix [A] 
In order to be as self-contained as possible, in Section [3] we present the group-theoretical 
backdrop and leave for the Appendix |B] a succinct exposition of the CWT on a general 
manifold X, collecting the main definitions used in this paper. In Section H] we briefly 
remind the CWT on K and extend it to the Lobachevsky plane C + and the open unit 
disk Di. The action of the affine group on C + extends naturally to the conformal group 
50(1,2) of the time axis KL This will serve us to introduce and establish a parallelism 
between standard and conformal wavelets in complex Minkowski space in Section We 
shall eventually work in the Cartan domain D 4 , although we shall provide in Section 15731 
an (intertwiner) isometry between the Hilbert spaces of holomorphic functions on O4 and 
the future tube domain C^_. The A-extended SMT turns out to be a valuable mathemat- 
ical tool inside D4 for proving admissibility and tight-frame conditions, reconstruction 
formulas and reproducing (Bergman) kernels in Section 15.21 We also discuss symmetry 
properties and comment on the Euclidean limit of the proposed wavelets in Sections 15.41 
and 15.51 respectively. Section |5J is devoted to convergence considerations and Section [7] to 
conclusions and outlook. In Appendix [C] we prove orthonormality properties of a basis of 
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homogeneous polynomials introduced in Section |5j 



2 Schwinger's Master Theorem: an Extension 

Schwinger's inner product formula [19] can be stated as follows: 

Theorem 2.1. (SMT). Let X be any 2x2 matrix X and Y = tl , where t is an arbitrary 
parameter and I stands for the 2x2 identity matrix. Let us denote by 



v j ( v~\ / (i + gQKi -gi) ! >p (j + ( j — q 2 



min(j+<ji,j+<j 2 ) 

( J + <l -\( 

'■ uu'xi 0,51+92) x " ' x ^ ^ 

v ™fc „j+m-fe„j+92-fc k- qi -q 2 fn\ 



the Wigner's T> -matrices for SU(2) JEE/, where j G N/2 (the spin) runs on all non- 
negative half-integers and qi,q 2 = —j, — j + 1, . . . , j — 1, j. Then the following identity 
holds: 



e (d u :X:d v ) e (u:Y T :v) 



= E ^" £ Z&W = det(I - iX)- 1 (4) 

j6N/2 g=-j 

where we denote by {u : X : v) = uXv T = £\ . =1 UiXijVj and d Ui = d/dui. 

This formula turns out to be essentially equivalent to MacMahon's Master Theorem: 

Theorem 2.2. (MMT) Lei X be an N x N matrix of indeterminates Xij, and Y be the 
diagonal matrix Y = diag(yi, y 2 , ■ ■ ■ , Vn)- Then the coefficient of y a = y^y^ 2 ■ ■ ■ y% N in 
the expansion o/det(7 — XY)~ l equals the coefficient of y a in the product 

N 

U(z<i2/i + x i2 y 2 H h x iN y N ) a ' . (5) 

8=1 

These abovementioned coefficients can be written in terms of the so-called SU(N) solid 
harmonics V^JX) (see [2U [25] and Appendix lAl for a general definition). SU(N) solid 
harmonics (11031) are a natural generalization of the standard Wigner's ©-matrices ([3]) 
to matrices X of size N > 2. In fact, replacing tX with XY in (J3J) and using the 
multiplication property 

E 2^w^(y) = ^(AT) (6) 

q'=-j 

and the transpositional symmetry 

ka y )= v3 ^ yT )i ( 7 ) 
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we can restate MMT for N = 2 as: 

E E ^W^(n = det(/-XF)- 1 . (8) 

jeN/2g,g'=-j 

Actually, MMT preceded Schwinger's result by many years. Schwinger re-discovered the 
MMT in the context of his generating function approach to the angular momentum theory 
of many-particle systems. The unification into a single form by using properties of the 
SU(N) solid harmonics was established by Louck in [23 1 12H 125] . 

Wigner matrices D 3 ,{X) are homogeneous polynomials of degree 2j in x^i- Inspired 
by Euler's theorem, we shall define the following differential operator: 

D x f(t) = (X + t^)f(t), AeN (9) 

which will be useful in the sequel. Now we are in condition to state and prove an extension 
of SMT 12.11 For the sake of completeness, a generalization for matrices X of size N > 2 
is also given in Appendix lAl 

Theorem 2.3. (A-Extended SMT) For every A G N, A > 2 and every 2x2 matrix X, 
the following identity holds: 

jeN/2 n=0 ^ ' ^ ' q=-j 

= det(J - tX)-\ (10) 

Proof: We start from the basic SMT 12.11 and apply the operator Di on both sides of 
Eq. flU): 

jeN/2 g=-j v y 

Here we have used that 

det(I - tX) = 1 - tr(tX) + det(tX) 

and that tr(Jf) and det(X) are homogeneous polynomials of degree 1 and 2, respectively. 
Making use of the expansion: 

J]t 2n det(X) n , (12) 



the expression ffTTj) can be recast as: 



£ (2j + 1) ^^ 2 «det(X)« £ ^ 9 (A) = ■ (13) 

jGN/2 n=0 q=-j ^ ' 



This identity is a particular case of (flOj) for A = 2. Now we shall proceed by induction 
on A. Assuming that ( TTOj) is valid for every A > 2 and applying the operator D\ on both 
sides of the equality ( flOjl . we arrive at: 

jeN/2 n=o ^ / V / 

xdet(xrt^w=A ; e -; 2 _^w , 

g=-i 



where we have made use again of (fl2l) . Considering ([121) one more time, we can assemble 
the previous expression as: 

E f^r E (" 1 - 2 2 ) (" + a - 2 ~ (A + 2J + 2")' 2)+2, " +m) 



jeN/2 n,m=0 



X ^X)T"^VUX)= ^ (14) 



i=-o 

Rearranging series: 

oo oo / n \ 

n,m=0 n=0 \m=0 / 

the identity (JT3J) can be recast in the form: 

jeN/2 n =0m=0 ^ ' ^ ' 

A 



x 

g=-j 

It remains to prove the following combinatorial identity: 



detW-g^W- w . (16) 



1 A An + A - 2\ /m +p + A - IV . 



A . 

m=0 



n + A — l\/n + p+ A 
A-l A A-l 



(17) 



We shall proceed by induction on n. Let us define both sides of the previous equality as 
the two sequences: 



„, . 1 (m + A - 2\ fm + p+ A - 1\ 



m=0 

n + A — 1\ /n+p + A 



A-l 7 V A-l 



It is easy to verify that F(0) = G(0). Assuming that F(n) = G(n), we ask whether 
F(n + 1) = G(n + 1). Indeed, on the one hand 



m=0 x ' x ' 

„, x A+p + 2(n + l) / n + A-l\ fn + p+X 
= b [nj H 



A — 1 V A — 2 7V A — 2 

v ; (A-l)(n+l)(n+p + 2) v ; 
(n + + p + 2) + (A + p + 2n + 2)(A - 1) 

(n + l)(n + p + 2) W 

and on the other hand 

v y \A-iyV A ~! / (n + l)(n + p + 2) v; 
Realizing that 

(n + l)(n + p + 2) + (A + p + 2n + 2)(A - 1) = (n + p + A + l)(n + A) 

we arrive at F(n + 1) = G{n + 1), which proves (iT7|) . Finally, inserting ( jTTl) in ( |T6i) . we 
conclude that ( |T0i) is valid for A + 1, thus completing the proof ■ 



3 The Group-Theoretical Backdrop 

The usual CWT on the real line R is derived from the natural unitary representation of the 
affine or similitude group G = SIM(1) in the space of finite energy signals L 2 (M, dx) (see 
Section H] for a reminder). The same scheme applies to the CWT on a general manifold X, 
subject to the transitive action, x — > gx, g G G, x G X, of some group of transformations 
G which contains dilations and motions on X. We address the reader to the References 
|H [5] for a nice and thorough exposition on this subject with multiple examples. For the 
sake of self-containedness, we also collect in the Appendix [B] some basic definitions which 
are essential for our construction of conformal wavelets. 

As already said in the Introduction, the CWT on spheres X = S N ~ l has been con- 
structed in [HI [9], [10] by means of an appropriate unitary representation of the Lorentz 
group in iV + 1 (space-time) dimensions G = SO(N, 1). The case of G = £0(2, 1) is 
particularly interesting as it encompasses wavelets on the circle S 1 and on the real line R, 
associated to the continuous and discrete series representations, respectively (see [27J for a 
unified group-theoretical treatment of both type of wavelets inside SL{2, R) ~ SO(l, 2)). 
The group SO '(1,2) (the conformal group in + 1-dimensions) has also been used to 
construct wavelets on the upper sheet of the two-sheeted hyperboloid H 2 [TT], or its 
stereographical projection onto the open unit disk (0Q). 
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The (angle-preserving) conformal group in N (space-time) dimensions is finite-dimensional 
except for N = 2. For N ^ 2, the conformal group SO(N, 2) consists of Poincare [space- 
time translations 6 M £ K. and restricted Lorentz £ SO + (N — 1, 1)] transformations 
augmented by dilations (a £ R+) and relativistic uniform accelerations (special confor- 
mal transformations c M £ R ) which, in A^-dimensional Minkowski spacetime, have the 
following realization: 

x'v = x >* + ft**, x'^ = A£(w)a^, 

_ ™^ _ Z M +C M Z 2 (18) 

X ^ ' X l + 2ciE + C 2 IE 2 ' 



iV-1 



respectively. We are using the Minkowski metric if = diag(l, — 1,^^, —1) to rise and 
lower space-time indices and the Einstein summation convention cx = c M x M . The new in- 
gredients with regard to the affine group SIM(1) are the extension from time-translations 
by 6° to A^-translations by & M , the addition of Lorentz transformations (rotations and 
boosts) and accelerations by & . Special conformal transformations can be seen as a 
sequence of inversions and translations by c M of the form: 

„ -4- T 2 nV 1 • ( o™/^ _1_ T ' / T 2 T*/^ -1- r^T 2 

~^ x 2 > x 2 y (x^ + x 2 c^) 2 /x 4 ~ 1 + 2cx + C 2 X 2 ' ^ ' 

They can also be interpreted as point- dependent (generalized/gauge) dilations in the 
sense that, while standard dilations change the spacetime interval ds 2 = dx^dxp globally 
as ds 2 — > a 2 ds 2 , special conformal transformations scale the spacetime interval point-to- 
point as ds 2 — > a(x)~ 2 ds 2 , with a(x) = 1 + 2cx + c 2 x 2 . The same happens with the 
squared mass m 2 , thus forcing a continuous mass spectrum unless m = 0, as for the 
electromagnetic field. 

The infinitesimal generators of the transformations f lT8|) are easily deduced: 



p — _JL ja — r _d _8_ 

n T n_d_ — — ?T T V -Q— 4- T 2 -Q— 

^ ~ X dxV A 4 9xV r X dxfi , 

and they close into the conformal Lie algebra 



(20) 



[M pi/ , M pa ] = r] up M^ a + f] Pi(J M V p - r) w M va - r] ua Mp P , 

[Pp,M pa \ = n w p a - v,*Pp, [p»p v ] = o, 

[K„,M pa ] = ri^Ka-V^Kp, [K^,K V ] = Q, (21) 

[D,Pp] = -P„, [D,Kp]=Kp, [D,Mp U ] = 0, 

[K„P V \ = 2{r 1 p V D + M llv ). 

The quadratic Casimir operator: 

C 2 = D 2 - X -Mp V M» v + \{PpK» + K^P% (22) 

generalizes the Poincare Casimir P 2 = PpP^ L (the squared rest mass). 

Any group element g £ 50(4,2) (near the identity element) could be written as the 
exponential map 

g = exp(ii), u = tD + b^Pp + c^Kp + w^M^, (23) 
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of the Lie- algebra element u (see J28J 122] )• The compactified Minkowski space M = 
§ N ~ l x Z2 S 1 , can be obtained as the coset M = SO(N, 2)/W, where W denotes the Weyl 
subgroup generated by K^, and D (i.e., a Poincare subgroup P = SO(N — 1, l)®!* 
augmented by dilations IR + ). The Weyl group W is the stability subgroup (the little group 
in physical usage) of x M = 0. 

For N = 2, the group 50(2, 2) is isomorphic to the direct product 5*0(1, 2) x 50(1, 2). 
It is well known that, in two dimensions, the conformal group is infinite dimensional. 
Actually, the splitting 50(2, 2) = 50(1, 2) x 50(1, 2) has to do with the separation into 
holomorphic and anti-holomorphic self-maps of the infinitesimal conformal isometries of 
a complex domain, the generators of which, 

L n = -z n+1 — , L n = -z n+1 — , z = x 1 + ix°, z = x 1 - ix°, neZ, (24) 
oz oz 

close into the Witt algebra [L m ,L n ] — (m — n)L m+n (idem for L). The conformal group 
in N = 2 + 1 dimensions, 50(3, 2), is also the symmetry group of the anti de Sitter space 
in 3 + 1-dimensions, AdS4 = 50(3, 2)/50(3, 1), a maximally symmetric Lorentzian mani- 
fold with constant negative scalar curvature (i.e., the Lorentzian analogue of 4-dimensional 
hyperbolic space) which arises, for instance, as a vacuum solution of Einstein's General 
Relativity field equations with a negative (attractive) cosmological constant (correspond- 
ing to a negative vacuum energy density and positive pressure). 

We shall focus on the 15-parameter conformal group in 3 + 1-dimensions, 50(4,2), 
which turns out to be locally isomorphic to the pseudo-unitary group 

SU(2, 2) = {ge Mat 4X 4(C) : g*Tg = V, det(<?) = 1} (25) 

of complex special 4x4 matrices g leaving invariant the 4x4 hermitian form V of sig- 
nature (+ H ). Here g^ stands for adjoint (or conjugate/hermitian transpose) of g (it 

is also customary to denote it by g*). Actually, the conformal Lie algebra (I2T1) can be 
also realized in terms of the Lie algebra generators of the fundamental representation of 
SU(2, 2), given by the following 4x4 matrices 



n mv- v = ^n"] _ i 

^ 2 1 4 4 



a u a^ 



= = I \ K » = 7 m1^7_ 5 / 



where 



2 V ) ' ' 2 V &» 



(26) 



r = I ^ I, 7 =*7 77 7 = I o 



denote the Dirac gamma matrices in the Weyl basis and 

"*-{\ ;)■•*-(! o)-*-(i ™ 
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are the Pauli matrices (we are writing = <x M ). Indeed, using standard properties of 
gamma and Pauli matrices, one can easily check that the choice (126]) fulfils the commu- 
tation relations ( J2TT) . 

To be more precise, SU(2,2) is the four-cover of 5*0(4,2), much in the same way as 
SU(2) is the two-cover of 50(3). This local isomorphism between the conformal group 
SO(N, 2) and the pseudo- unitary group SU(M, M) only happens for N = 1 and iV = 4 
dimensions, where 

50(1, 2) = SU(1, 1)/Z 2 , 50(4, 2) = SU(2, 2)/Z 4 . (28) 

The A-extension of the SMT given in Theorem 12.31 (and Theorem IA.2j) turns out to 
be closely related to the group SU(2, 2) (and SU(N,N) in general), providing a kind 
of generating function for the unitary-representation functions of this group (the discrete 
series, to be more precise). This formula will be a useful mathematical tool for us, specially 
in proving admissibility and tight-frame conditions and providing reconstruction formulas. 
From this point of view, the conformal group SO(N,2) in iV = 4 dimensions is singled 
out from the general iV-dimensional case, at least in this article. 

Before tackling the construction of conformal wavelets in 8-dimensional complex Minkowski 
space in Section [5J we shall briefly remind the simpler case of the CWT on the time axis 
R and its extension to the Lobachevsky plane C + and the open unit disk © 1; which are 
homogeneous spaces of 5*0(1, 2). 

4 Wavelets for the Affine Group 

Let us consider the affine or similitude group of translations and dilations in one dimen- 
sion, 

G = SI M{\) = RxR + = {g = (b,a)/beR,ae R + }, 
with group law (g" = g'g): 

a" = a' a 
b" = b' + a'b 

This group will serve us as an introduction for studying the most interesting case of 
the conformal group G = 50(4,2) as a "similitude" group of space-time, which will be 
considered in the next section. 

The left-invariant Haar measure is: 

du(g) = —da A db 
a 2 

The representation 

[U x (a,b)</>](x) = a~ A 0(^— -) = <Mx) (29) 

a 

of G on L 2 (R, dx) is unitary for A = ~ + is. In fact, every U\ is unitarily equivalent to 
Ux/2 and one always works with A = |. This representation is reducible and splits into 
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two irreducible components: the positive u > and negative u < frequency subespaces. 
Restricting oneself to the subspace u > 0, the admissibility condition f |123[) assumes the 
form 

I du < ^ 

Jo w 

where ^ stands for the Fourier transform of if). Given an admissible function if) G 
L 2 (M, dx), the machinery of wavelet analysis proceeds in the usual way. 

4.1 Wavelets on the Lobachevsky plane C + 

An extension of the representation ( 1291) of the affine group, this time on the space 
L\(<C + , di)x) of square integrable holomorphic functions on the upper half complex plane 
(or forward tube domain) 

Ti = C+ = {w = x + iy G C/Q(w) = y > 0}, (30) 

is given by: 

[U x (a,b)<f>)(w) = a- x <f>C^). (31) 

a 

This representation is unitary with respect to the scalar product: 

r 2A — 1 

W)= <f)(w)<j>'(w)dv x (w,w), du x (w,w) = — Qt(w) 2 W\dw\, (32) 

Jc+ 47r 

for any 0,0' G L\(£.j r ,dP\), where we use \dw\ as a shorthand for the Lebesgue measure 
d$l(w) A d$$(w). Although all representations U\, A > 1, are equivalent, they become 
inequivalent when the affine group is immersed inside the conformal group of the time 
axis R, SO(l, 2) ~ SL(2, R) ~ SU(1, 1). Actually, this will be the case with the conformal 
group 50(4, 2) in the next section. This immersion of SIM(1) inside SL(2, R) is apparent 
for the Iwasawa decomposition KAN (see, for instance, [30J) when parameterizing and 
element g G SL(2,R) as: 

a= (co S 6 -sin6\(j= a \fl b\ = (<% ^-v^sin£\ 
y \ sin9 cos9 J \ ^ J \ 1 J v4cos^ + ^ J 

where 9 G (— 7r, tt] (see [27] for a unified group-theoretical treatment of wavelets on M and 
the circle S 1 inside SL(2,R)). 

4.2 Wavelets on the open unit disk Di 

There is a one-to-one mapping between the Lobachevsky plane C + and the open unit disk 
(or Cartan domain) 

©i = {z G C, \z\ < 1}, (33) 
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given through the Cayley transformation: 

. . 1 + iw . s 1 — z 

z(w) = <H- w(z) = i . (34) 

v ' l-vw v ' 1 + z v/ 

Note that the (Shilov) boundary S 1 = {z G C : \z\ = 1} of Di is stereographically 
projected onto the boundary K = {w G C : ^s(w) = 0} of C + by w(e 10 ) = tan(#/2). 

We can establish an isometry between L\(C + , du x ) and the space L\(JD>i, dis\) of square 
integrable holomorphic functions on the unit disk Di with integration measure 

2A — 1 

du x (z,z) = — (1 - zzf^\dz\, A > 1, (35) 

7T 

where z denotes complex conjugate. This isometry is given by the correspondence 

S x : L 2 h {B u du x ) — ► Ll(C+,dv x ) 
cj) i — > S X (f) = <j>, 

with 

4>{w) = 2 2A (1 - iw)- 2X (p(z(w)) (36) 

and z(w) given by (l34"j) . In fact, taking into account that (1 — zz) = 2 2 Q(w)\l — iw\~ 2 
and the Jacobian determinant |g?z|/|g?k;| = 2 2 |1 — iw\~ A , then 

o \ _ I r 

%|(b = / <t>{z)<t>'{^-^f {X - l) \dz\ 

2A - 1 



47T 



2 2A(! _ t W )-^ z ( w ))2 2 \l - iwy^^'iziw^iwf^ldwl 
<p(w)<p(w)dv x (w,w) = (0|0 / ) L 2(c+)- 
The constant factor (2A — l)/7r of du x (z, z) is chosen so that the set of functions 

1 /2 

Vn {z)= ( 2X + n ~ 1 ) z n y n = ,l,2,..., (37) 



n 

constitutes an orthonormal basis of L^(Bi, dv x ), as can be easily checked by direct com- 
putation. These basis functions verify the following closure relation: 



J2< Pn (z)<p n {2f) = (l-Z2f)- 2X , (3? 



n=0 



which is nothing other than the reproducing (Bergman) kernel of L\(Jix,dv\) (see e.g. 
[4 J for a general discussion on reproducing kernels). We shall provide a four- dimensional 
analogue of (157)1 and ( 155|) in Eqs. (155]) and (jSTJ, respectively, and prove the orthonormality 
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in Appendix O The isometry S\ given by (1361) maps the orthonormal basis (1371) of 
L 2 h (Pi,dv\) onto the orthonormal basis 

<p n (w) = 2 2A (1 - m)- 2 V(^(^)), n = 0, 1, 2, . . . (39) 
of L^C.)-, dz/), which verify the reproducing kernel relation 

oo 

^^Hy5„K) = C-(w - w'))- 2 \ (40) 

n=0 

Let us denote by Va = S^UxSx the representation of the affine group on L\(J5)\,dv\) 
induced from f[3~Tj) through the isometry f[3"oT) . More explicitly: 

[V A (a,6)0]^) = [S?Ux(a,b)<j>)(z) = a' A 1 — -±- J f * f 
This representation is, by construction, unitary on L^(Di,d^). 



5 Wavelets for the Conformal Group 50(4, 2) 

The four-dimensional analogue of the extension of the time axis K to the time-energy 
half-plane C + is the extension of the Minkowski space M 4 to the (eight-dimensional) 
future tube domain C 4 of the complex Minkowski space C 4 (see later on this Section). 
The four-dimensional analogue of the one-to-one mapping between the half-plane C + 
and the disk Di is now the Cayley transform fj471) between C 4 and the Cartan domain 
B 4 = U(2,2)/U(2) 2 , the Shilov boundary of which is the compactified Minkowski space 
£7(2) (the four- dimensional analogue of the boundary U(l) = S 1 of the disk Dx). Let us 
see all this mappings and constructions in more detail. 

5.1 Wavelets on the forward tube domain 

The four- dimensional analogue of the upper-half complex plane ( 130]) is the future/forward 
tube domain 

T 4 = C 4 = {W = X + %Y = G Mat 2x2 (C) : Y > 0} (41) 

of the complex Minkowski space C 4 , with X = x^cr^ and Y = y^a^ hermitian matrices 
fulfilling the positivity condition Y > <^ y° = Q(w°) > \\y\\. 

The domain C 4 is naturally homeomorphic to the quotient SU(2,2)/S(U(2) x £7(2)) 
in the realization of the conformal group in terms of 4 x 4 complex (block) matrices / 
fulfilling 

G = 5/7(2,2) = if=(* T l ^\e Mat 4x4 (C) : fTf = I\det(/) = l} , (42) 



13 



with 

I 

1 



7° 



the time component of the Dirac 4x4 matrices 7 M in the Weyl basis (/ = cr° is the 2x2 
identity matrix and /' = /* stands again for adjoint or conjugate/hermitian transpose of 

/). In general, T is a 4 x 4 hermitian form of signature (+ H ). The inverse element 

of / is then given by: 

/-■=tw=($ if 

The particular identification of with the coset SU(2, 2) / S{U{2) 2 ) is given through: 

W = W(f) = (S + iR)(Q + iT)- 1 = (Q + iTy^S + iR). (43) 

The left translation /' — > ff of G on itself induces a natural left-action of G on C 4 given 
by: 

W = W(f') 4f' = W(ff') = (RW + S){TW + Q)-\ (44) 
Let us make use of the standard identification x M -<-»• X = x M cr M between the Minkowski 



space M 4 and the space of 2 x 2 hermitian matrices X, with er M the Pauli matrices (J2j 
and M = det(X) the Minkowski squared-norm. Setting W = the transfor- 

mations (118j) can be formally recovered from (1441) as follows: 

i) Standard Lorentz transformations, x' M = A$(u))x u , correspond to T = S = and 
R = Q~ 1Jf G SL(2, C), where we are making use of the homomorphism (spinor map) 
between S0+(3, 1) and SL(2, C) and writing W = RWR\ R G SL(2, C) instead of 
x'» = A£x v . 

ii) Dilations correspond to T = S = and R = Q~ x = a 1 / 2 / 

iii) Spacetime translations correspond to R = Q = I and S = b^a^, T = 0. 

iv) Special conformal transformations correspond to R = Q = I and T = c M <r M , S = 
by noting that det(J + TW) = l + 2cx + c 2 x 2 . 

We shall give the next proposition without proof. Instead, we address the reader to 
its counterpart (Proposition 15.21) in the next Subsection for an equivalent proof. 

Proposition 5.1. The representation of G on square-integrable holomorphic functions 
f(W) given by 

[U x {f)(f]{W) = detGRt - T*W)- X (p{{Q^W - S^)(R^ - T^W)' 1 ) (45) 
is unitary with respect to the integration measure 

dP x (W,W^) = ^det(i(W*- W)) x - 4 \dW\ = ^(w) 2 ^\dW\, (46) 
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where A G N, A > 3 (the "scale or conformal dimension") is a parameter labelling non- 
equivalent representations, c\ = (A — 1)(A — 2) 2 (A — 3)/7r 4 and we are using \dW\ = 
/V i j- Q dM(w ti )dSl(Wf t ) as a shorthand for the Lebesgue measure on C+. 

We identify the factor M{f,Wf/ 2 = det( J R t - T^W)~ X in @5) as a multiplier or 
Radon-Nikodym derivative, (remember the general definition in ( 1121D ) . It generalizes the 
factor a~ x in f l3Tj) by extending (global) standard dilations i? = Q^ 1 = a 1 / 2 /, T = 5 = 
to (local/point-dependent) "generalized dilations" with T = c^o^ . The representation 
f )45|) is a special (spin-less or scalar) case of the discrete series representations of SU(2, 2), 
which are characterized by A and two spin labels s\ and s 2 . Decomposing the discrete 
series representations of SU(2, 2) into irreducible representations of the inhomogeneous 
Lorentz group leads to a continuous (Poincare) mass spectrum [T7j. 

5.2 Wavelets on the Cartan domain D4 

Instead of working in the forward tube domain C+, we shall choose for convenience a 
different eight-dimensional space O4 generalizing the (two-dimensional) open unit disk 
Di in fl33l) . where we shall take advantage of the full power of the A-extension of the 
SMT given by the formula ffTUj) . Both spaces, and ©4, are related by a Cayley-type 
transformation, which induces an isomorphism between the corresponding Hilbert spaces 
of square-integrable holomorphic functions on both manifolds (see later on Section I5T31 . 

5.2.1 Cayley transform and D 4 as a coset of S77(2,2) 

The four-dimensional analogue of the map (1511) form the Lobachevsky plane C + onto the 
unit disk Di is now the Cayley transformation (and its inverse): 

W -> Z(W) = (I-iW)- 1 (I + iW) = (I + iW)(I-iW)- x , , * 

z -> w(z) = i{i - z){i + zy 1 = i(I + Z)-\l-Z), { 7) 

that maps (one-to-one) the forward tube domain C 4 onto the Cartan complex domain 
defined by the positive-definite matrix condition: 

D 4 = {Z e Mat 2x2 (C) : I - ZZ^ > 0}. (48) 

Note that the (Shilov) boundary 

D 4 = 1/(2) = {Z e Mat 2x2 (C) : ZZ ] = Z f Z = 1} = § 3 x Za S 1 

of D4 is a compactification of the real Minkowski space 

M 4 = {W G Mat 2x2 (C) : = W} } 

i.e., the boundary of C 4 (see e.g. [IB])- The restriction of the Cayley map Z — > W(Z) to 
Z G U(2) is precisely the stereographic projection of U{2) onto M 4 . 
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The Cartan domain O4 is naturally homeomorphic to the quotient SU(2, 2)/S(U(2) 2 ) 
in the new realization of: 

G = SU(2,2) = {0 = ( £ J e Mat 4x4 (C) : g^g = ^,det(g) = lj , (49) 



with 

,5 / -/ 



7 - -1 / 



the fifth Dirac 4x4 gamma matrix in the Weyl basis (J = er° denotes again the 2x2 
identity matrix). The inverse element of g is now: 



9 1 = tVt 5 



^ -ct 

-fit £>t 

Both realizations f l42|) and f l49|) are related by the map 

f -( A B \ - - l ( R + 1S-1T + Q -R + 1S + 1T + Q , .. 

! ^ g ~\C D)~ l 2 \ —R - iS - iT + Q R-1S + 1T + Q ' ' ( ' ,U) 

with 

t = -LC 7 -' 

The particular identification of ©4 with the coset SU(2, 2) / S{U{2) 2 ) is given by [see later 
on Eq. (}59~]) for more details]: 

Z{ J g) = BD-\ Z\g)=CA~\ (51) 

Actually, making explicit the matrix restrictions g^ 5 g = 7 5 in (149]) : 



£>tD - fit 5 = / 
9-^ = J 4 x4 ^ <( ^ A - CtC = / (52) 
A^B-CW = 

and 

r £>£>t - CCt = / 

gg- l = h x4 ^\ AA^-BB^ = I (53) 

[ A& - BD^ = 0, 

the positive matrix condition in (148]) now reads 

I-Ztf = I- A- l] C ] CA~ l = (AA^y 1 > 0, (54) 

where we have used the second condition in (152]) . Moreover, using the identification (|5ip 
and the first condition in (152]) . we can see that 

det(ZZt) = det(fitfi) det(J + fitfi) -1 < 1, (55) 
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This determinant restriction can also be proved as a direct consequence of the positive- 
definite matrix condition / — ZZ^ > 0. In fact, the characteristic polynomial 



det((l- p)I - ZZ r ) = 1 -%x{pI + ZZ^) + det(p/ + ZZ t ) 

= p 2 -(2-tr(ZZ t ))p + det(J -ZZ ] ) 

yields the eigenvalues 

P± = 2 ~ tr( ^ t)±V ^ , A = (2 - tr(ZZt)) 2 - 4det(/ - Z&). (56) 

Since / — ZZ* is hermitian and positive definite, its eigenvalues p± are real and positive. 
The condition p_ > implies that: 

2-ti{ZZ ] ) > ^tr(ZZ^) < 2, (57) 

and the fact that A > gives: 

< A = ti(ZZ^) 2 -Adet(ZZ^) &et(ZZ ] ) < hr(ZZ^) 2 < 1, (58) 

where we have used (1571) in the last inequality. From (1571) . we can regard O4 as an open 
subset of the eight-dimensional ball with radius Vz. All those bounds for ZeD 4 will be 
useful for proving convergence conditions later on Section [6j See also the Appendix ICl for 
a suitable parametrization of Z when computing scalar products. 

5.2.2 Haar measure, unitary representation and reproducing kernel 

Any element g e G admits a Iwasawa decomposition of the form 

( A B \ ( A 1 ZA 2 \(U 1 0\ 
9 = ( C D J = ( ZtA a A 2 J ( £7 r ) ' (59) 



2 

with 

A 1 = (AA ] ) l/2 = (I- ZZ ] y l/2 , Ui = A' 1 A, 
A 2 = (DD ] ) l/2 — (I — Z^Z)- 1/2 , U 2 = A 2 l D. 

This decomposition is adapted to the quotient D 4 = G/H of G = SU(2, 2) by the maximal 
compact subgroup H = S(£7(2) 2 ); that is, Ui,U 2 G £7(2) with det(UiU 2 ) = 1. In order 
to release U\ y2 from the last determinant condition, we shall work from now on with 
G = £7(2,2) and H = £7(2) 2 instead. Likewise, a parametrization of any £7 e £7(2), 
adapted to the quotient S 2 = £7(2)/£7(l) 2 , is (the Hopf fibration) 
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where z = b/d G C ~ § 2 (the one-point compactification of C by inverse stereographic 
projection), 5 = (1 + zz)- 1 ' 2 and e i/Sl = a/|a|, e ifti = d/\d\. 

The left-invariant Haar measure [the exterior product of left-invariant one-forms g~ l dg\ 
of G proves to be: 

dfi(g) = d/i(g)\ G/H d/i(g)\ H , (61) 
dfi(g)\ G/H = det(I-ZZ^)- 4 \dZ\, 
dfi{g)\ H = dv{U 1 )dv(U 2 ), 

where we are denoting by dv(U) the Haar measure on £7(2), which [using (I60j) ] can be in 
turn decomposed as: 

dv(U) = dv(U)\ u(2yu{1)2 dv{U)\ u{1)2 , 

MU)\ U(2)/U{1) 2 = dv(U)\ s2 = ds(U) = (l + zz)- 2 \dzl (62) 
dv(U)\ u(l)2 = dfadfo, 

where \dz\ and \dZ\ denote the Lebesgue measures in C and C 4 , respectively. 
Let us consider the space of holomorphic functions <f>(Z) on D4. 

Proposition 5.2. For any group element g G G, the following (left-)action 

cj> g (Z) = [U x (g)<p](Z) = fet(D* - Z)~ X <I>(Z'), (63) 
Z' = g- x Z = (A^Z -C ] )(rf -B^Z)- 1 

defines a unitary irreducible square integrable representation of G on L^(B4, dv\) under 
the invariant scalar product 



>>') = / c(>(Z)cl>'(Z)du x (Z,Z^, 

du x (Z,Z^) ee c x det(I - ZZ^) x - 4 \dZ\, (64) 

for any A G N, A > 4 (the "scale or conformal dimension"), where c\ = 7r _4 (A — 1)(A — 
2) 2 (A — 3) is chosen so that the unit function, 4>(Z) = 1, WZ G ©4, is normalized, i.e. 

(010) = 1- 

Proof: One can easily check by elementary algebra that U\(g)U\(g') = U\(gg'). In order 
to prove unitarity, i.e. (0 9 |0 S ) = (0|0) for every g G G, we shall make use of (1521) and 
(P|) . In fact: 

det(I-Z'Z^) = |det(D f - B^Z)\~ 2 det(I - ZZ^), 
and the Jacobian determinant 

\dZ\ = \dZ'\\det(D^ - B^Z)\ 8 , 

give the isometry relation ||0 9 || 2 = ||0|| 2 . Now taking g' = g~ x implies the unitarity ofU\. 
For the computation of c\ and other orthonormality properties see Appendix O B 
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In the next Section, we shall provide an isomorphism between L^D^, dv) and L|(Ct, dv\ 
where we enjoy more physical intuition. 

In order to prove admissibility conditions in Section 15.2.31 it will be convenient to give 
an orthonormal basis of L^(B4, dv\). 

Proposition 5.3. The set of homogeneous polynomials of degree 2j + 2m: 

< m J z ) - \fWn^) FSF) det(zrvuz), ( 65 ) 

m G N, j G N/2, qi , q 2 = -j, + 1, j, 

constitutes an orthonormal basis of L 2 h {0^dux), that is: 

(<r 92 I^S) = 5 i,i' 5 rn, m '5 qi , q ' i 5 q2A > 2 . (66) 

Note that the number of linearly independent polynomials JJ? J=1 z™- 3 of fixed degree 

of homogeneity n = Y2fj=i n ij * s ( n + 1) (w + 2) (n + 3)/ 6, which coincides with the number 
of linearly independent polynomials ( 1651) with degree of homogeneity n = 2m + 2j. This 
proves that the set of polynomials (165]) is a basis for analytic functions 4> G L^(B 4 , dz/ A ). 
Moreover, this basis turns out to be orthonormal. We address the interested reader to 
the Appendix [C] for a proof. We prefer to omit it here in order to make the presentation 
more dynamic. 

Note also the close resemblance between the definition ( 16"5]) and the left-hand side of 
the equality ( TTOj) in the A-extended SMT I2H1 In fact, taking tX = Z^Z' in (flOj) and using 
the properties © and (jTj) of Wigner's D-matrices, we can prove the following closure 
relation for the basis functions (1651) : 

oo j 

EEE ^(Z)^(Z') = det(7 - tfZ<)-\ (67) 

jSN/2 m=0 q,q'=-j 

which is nothing other than the reproducing (Bergman) kernel in L^(B 4 , di>\) (see e.g. [I] 
for a general discussion on reproducing kernels). Note that, although the scalar product 
(|64p is only valid for A > 4, the expression fIBTI) is formally valid for A > 2, since we are 
just using in it the requirements of the A-extended SMT 12.31 The case A = 2 is related to 
the Szego kernel (see e.g. [IS]). 

5.2.3 Admissibility condition, tight frame and reconstruction formula 

Theorem 5.4. The representation UJ3) is square integrable, the constant unit function 
ip(Z) = <£>q'q(Z) = 1, \/Z G ©4 being an admissible vector (fiducial state or mother 
wavelet), i.e.: 

c^= [ \(U x (g)^)\ 2 dfi(g) < oo (68) 
Jg 
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and the set of coherent states (or wavelets) F = {ip g = U\(g)ip,g G G} constituting a 
continuous tight frame in L^iJDi^, du\) satisfying the resolution of the identity: 



A= / \^ g ){^ g \dfi(g) = c^l. (69) 

J G 

Proof: Using the extended SMT [2~3l for tX = D~ X CZ, we can expand 

ip g (Z) = det( J D t - B ] Zy x = det{D j )- x det(J - (BD- l ) ] Z)- x 

j=0 n=0 v 7 v 7 

q=-j 

Now, taking into account that det((B D- l ) ] Z) n = det(( J B£)- 1 ) t ) n det(Z) n and the prop- 
erty ([6]) for 

Vl q ((BD-yZ) = £ Vl ql ((BD-y)Vl, q (Z), (71) 
q'=~j 

we recognize the orthonormal basis functions (l65i) in the expansion ( ITU]) , so that we can 
write the coherent states (wavelets) as: 

oo j 

uz) = E E E ( 72 ) 

jSN/2 n=0 q,q'=-j 

with "Fourier" coefficients 



X 

q=-j 



E ^V^ 1 )^ = det(D)-^ q r q (BD-±). (73) 



Using the orthogonality properties fl66l) of the basis functions fl65l) . we can easily compute 
|(^ A (^)^|^>| 2 = K^l^)! 2 = l<'o(^)| 2 = det(,DZ?t ) -A = det(/ _ ^ +) a ? (74) 

where we have defined Z = BD~ X and used the first condition in (|52p . Using the Haar 
measure ( 16T|) . the admissibility condition (168]) gives: 

c^= / c//i(^)| G/H det(/-ZZ t ) A / dv(U x )dv{U 2 ) = c^ 1 (^-] <oo, (75) 
Jg/h 
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where we have identified dfi(g)\ G ^ H det(7 — ZZ^) X = c^dv\{Z, Z^) and taken into account 
that J B dv\(Z, Z^) = 1 and 

(27r times the area of the 3-sphere § 3 = SU(2) of unit radius). 

Now it remains to prove that the resolution operator f[6"9"j) is a multiple of the identity 
I in L|(B 4 , dv\). For this purpose, we shall compute its matrix elements: 



= v{U{2)f I d^g)\ G/H det(I - ZZ^™ 2 (Z)^(Z) 

J G I H 

= c^5j ) ji5 m)m i5 qitq ' i 5 q2:q ' 2 , (77) 

where we have used (1731) . the orthogonality properties ( J66l) of the basis functions (!65l) and 
the fact that G/H = D 4 . ■ 

The reconstruction formula (11281) here adopts the following form: 



G* 



0(Z) = / $Ml> g (Z)dn(jg), (78) 

with wavelet coefficients 



<Ms) = -(^10) = — / det(D - &B)- x <p(Z)dv x (Z, tf). (79) 



Expanding in the basis ([65 



OO J 



w = E E E C'C'(^)> 

jSN/2 n=0 q,q'=-j 

and using (!72|) together with the orthogonality properties (166|) . we can write the wavelet 
coefficients (1T9"j) in terms of the Fourier coefficients 0^,™ as: 

<m*) = - E E E = - E E E ^(/^v^zr^-. 

^ jeN/2 m=0 q,g'=-j ^ jeN/2 m=0 q,q'=-j 

5.3 Isomorphism between L|(D4,(iz/) and L|(C^,rfz>) 

For completeness, we shall give an isometry between L^(©4,c?z/) and L\(C\,du) which 
allows us to translate mathematical properties and constructions from one space into the 
other. 
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Proposition 5.5. The correspondence 

S x : L 2 (3 4 ,du x ) — > L 2 (C%,dv x ) 

4> 1 — > <S\4> = 4>, 

with 

$(W) = 2 2X det(J - iWy x <p(Z(W)) (80) 
and Z(W) given by the Cayley transformation [J^D , is an isometry. Actually 

(0|0')l2(d 4 ,^ a ) = (S\<P\S\<P') Ll{c\,dvxY ( 81 ) 

Moreover, S x is an intertwiner (equivariant map) of the representations ( Tffgj) and fT^l ). 
that is: 

U x = S^U X S X . (82) 
Proof: The left-hand side of the equality flHTj) is explicitly written as: 

'W,^) = / W)<P\Z)c x det{I - Ztf) x - A \dZ\. 

Taking into account that 

det(J - ZZ ] ) = det(2z(W /t - W))\ det(/ - iW)\- 2 
and the Jacobian determinant 

\dZ\ = 2 12 \det{I -iW)\~ 8 \dW\, 

then 

dv x {Z,Z x ) = c A det(J- ZZ ] ) x ~ 4 \dZ\ = 

= 2 4A " 4 | det(J - tW)\~ 2X c x detfyw* - W)) x ~ 4 \dW\ 
= 2 4A |det(J -iW)\~ 2X dv x (W,W ] ), 

which results in: 



<f>(z)<f>(z)di, x (z,tf)= / 4>(w)4>(w)di> x (w,w^, 



C 4 



thus proving flHTl) . 

The intertwining relation ( 182]) can be explicitly written as: 



[W A 0](Z) = det(L>t - £t£)-A0((4t£_ct)(£)t _£t£)-i) 



3^ 



(Z) = det(J-iIV) A det( J R t -Ttiy)- A det(/-ziy , ) _A 0(^(W)) ; 



(83) 
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where W = (Q^W — S')(R* — T^W) 1 . On the one hand, we have that the argument of 



is: 



Z{W') = (I + zW')(I — iW')~ l 

= ((fit - T*W) + i(Q ] W - S ] )) ((fit - T*W) - i(Q ] W - 5 t )) _1 
= ((fit - i£fl) + i{Q ] + iT r ) W) ((fit + itf) - i{Q ] - iT*)W) _1 . 



-i 



Taking now into account the map (150]) we have: 

Z(W) = ((^ - &) + + Cft)W) ((fi f - fit) - ipt + fit)^)" 1 
= (At(/ + zW) - Ct(J - zW)) (fit(J _ iw) - B\I + zW)) 
= (A^Z-C^)(D^-B^zy\ 

as desired. On the other hand, we have that 

(/ - iW'){R) - T^W) = (fit _ T^W) - i{Q ] W - S ] ) = (fit + itf) - i{Q ] - iT r )W 
= (fit _ B r ) - i(fit + B ] )W = D\l - iW) - fit(J + iW) = (fit - B ] Z)(I - iW) 

which implies 

det(I - iW) x det(fit - fit^)- A det(J - iW'y X = det(fi f - fit^)~ A (84) 

That is, the equality of multipliers in (1831) . ■ 

As a direct consequence of Proposition I5.5[ the set of functions defined by 

« 2 W = 2 2X det(I -iW)- x ^ q2 (Z(W)), (85) 

with <^q™ q2 defined in (165|) . constitutes an orthonormal basis of fi|(C* , g?z>a) an d the closure 
relation 

oo j 

EEE ^r q (w)^r q (w') = detC-(w^-w')r\ m 

jeN/2 m=0 q,q'=-j 

gives the reproducing (Bergman) kernel in L|(C+, g?i>a)- 

The isometry ( )80|) also allows us to translate the results of Theorem l5.4l form L^(D 4 , dv\) 
into fi|(C+, c?i> A ). Indeed, from (18"0"|) we conclude that the function ^ G L\(C^, <iz>A) given 
by: 

^(jy) = 2 2A det(/-zW)- A (87) 

is admissible. The construction of a tight frame and a reconstruction formula form this 
mother wavelet parallels ( )69|) and ( !7H|) . respectively. 
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5.4 Symmetry properties of the proposed conformal wavelets 



When working with wavelets on the sphere [U |9j [10] it is customary to take axisymmetric 
(or zonal) wavelets, that is, admissible vectors ip which are invariant under rotations 
around the (namely) z-axis, although more general implementations including directional 
spherical wavelets are also possible (see e.g. [31]). Let us discuss the symmetry properties 
of our proposed admissible wavelets flH7|) . Applying a general SU(2, 2)-transformation 
(USD to (|87D gives: 

[U\(f)$\(W) = 2 2A det(i?t _ TWy x det(I - iW')-\ 

W' ={Q^W -St)(i?t _ T ] W yi V ) 

Using the identity ( 1841) we have: 

[Ux(f)i>](W) = det(,D t - BiZ)- x <$(W), (89) 
which leaves invariant ip (up to a global phase) if: 

B = 0^C = 0^S = -T, Q = R, (90) 
where we have used (|52|) . ( 153]) and ( 150]) . Thus, the elements / G SU(2, 2) of the form 

f=( H l ») (91) 



x iS R , 

leave invariant ( 1H7|) . The constraints /'7°/ = 7 imply: 

+ R)R = /, srti? = tfs. (92) 

For 5 = 0, R is unitary. For R = I, S is hermitian with S* 2 = 0. The last condition is 
satisfied for translations S = b^a^ along null (light-like) vectors b 2 = b^ = det(S') = 0. 
This leaves us a 7-dimensional subgroup of SU(2,2), isomorphic to S(U(2) x U(2)), as 
the isotropy subgroup of the admissible vector ( 1871) . Any other basis state ( |85l) could be 
used as a fiducial state to construct oriented wavelets. 

In the next Figure we provide a visualization of this wavelet (modulus and argument) 
for the particular case of W = wa° (temporal part), w = x + iy, for which ip(W) = 
2 2A (1 — iw)~ 2X reduces to <fo(w) in f )39|) . We take A = 1 for simplicity. 
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5.5 The Euclidean limit 



We have seen that the Shilov boundary of O4 is the compactified Minkowski space U(2) = 
§ 3 x% 2 S 1 (the four- dimensional analogue of the boundary U(l) = S 1 of the disk Di). One 
expects the wavelet transform on to behave locally (at short scales or large values of 
the radius p) like the usual (flat) wavelet transform on M. N . Indeed, in [27], one of the 
authors and collaborators discussed the Euclidean limit (infinite radius) for wavelets on 
S 1 . The procedure parallels that of Ref. [i9J for wavelets on § 2 . In these references, the 
Euclidean limit is formulated as a contraction at the level of group representations. Let 
us restrict ourselves, for the sake of simplicity, to the conformal group SO(l,2) in 1+0 
(temporal) dimensions. The realistic 1+3 dimensional case 5*0(4, 2), although technically 
more complicated, follows similar guidelines and will be left for future work. 

Let us denote simply by P = P and K = K the temporal components of P p and K p 
(the generators of spacetime translations and accelerations). The Lie algebra commutators 
of 5*0(1, 2) are [remember the general iV-dimensional case (121j) ]: 

[D, P] = -P, [D, K) = K, [K, P] = 2D. (93) 

A contraction Q' of the Lie algebra Q = so(l, 2) along sim(l) (generated by P and D) can 
be constructed through the one-parameter family of invertible linear mappings n p : M 3 — > 
M 3 , p G [1, 00) defined by: 

ir p (D) = D, n p (P) = P, tt p (K) = p^K, (94) 

such that the Lie bracket of Q' is: 

[X,Y]'= lim TrrVp^VH, (95) 

p—^oo r 

with [*, •] the Lie bracket fl9"3"l) of Q. The resulting Q 1 commutators are: 

[D, P] = -P, [D, K}' = K, [K, P]' = 0. (96) 

The contraction process is lifted to the corresponding Lie groups G' = M 2 x 1R + and G = 
50(1, 2) by considering the exponential mapping e^f . The idea is that the representation 
of G contract to the usual wavelet representation of the affme group SIM(1) in the 
following sense: 

Definition 5.6. Let G' be a contraction ofG, defined by the contraction map U p : G' — > G, 
and letW be a representation of G' in a Hilbert space W '. Let {U p },p G [l,oo) be a one- 
parameter family of representations of G on a Hilbert space H, p , and i p : H, p — )■ T> p a linear 
injective map from 1-L P onto a dense subspace T> p C H' . Then we shall say that W is a 
contraction of the family {U p } if there exists a dense subspace V C TV such that, for all 
<fi G V and g' G G' , one has: 

• For every p large enough, G T> p and Up{U. p {g'))i~ l 4> G L p lr L) p . 
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. lim^ oo ||i p W p (n p (^))t-V-WV)0ll«' = O ) Wee. 
More precisely, one can prove that: 

Theorem 5.7. The representation [W(b,a)(j)](x) = ^</K~0 of the affine group SIM{1) 
is a contraction of the one-parameter family U p of representations of SO(l, 2) on "H p as 
p — > oo. That is: 

]M\\tMJn^(b,a))L- 1 <f)-U\b,a)(f>\\ w = 0, V(M)eR^R + , (97) 

where : SIM(1) 5*0(1, 2)/R is the restricted contraction map, with a' : G'/M. — > G' 
a given section. 

This construction can be straightforwardly extended to G — SO (4, 2), the contraction 
G' being the so-called G15 group of Ref. [32J. A thorough discussion of the Euclidean 
limit of the conformal wavelets constructed in this paper falls beyond the scope of this 
article and will be left for future work |33j. Here we just wanted to give a flavor of it. 



6 Convergence Remarks 

Schwinger's Theorem 12. II and its extension 12 .31 have been stated in the sense of generating 
functions in terms of formal power series in some indeterminates. From this point of 
view, we have disregarded convergence issues. However, infinite series expansions like for 
instance ( 112]) would require in particular that |t| 2 | det(X)| < 1. We shall prove that such 
convergence requirements, together with additional restrictions coming from the basic 
Theorem 12.11 are automatically fulfilled inside the complex domain O4 for tX = Z^Z, 
with Z = BD' 1 in the expansion flTOj) . Let us state these convergence requisites. 

Proposition 6.1. A sufficient condition for the convergence of the expansions Q) and 
(T72P for t = 1 is that: 

\x n \ < 1, \x 22 \ < 1, \xi 2 x 2 i\ < 1, I det(X)| < 1. (98) 



Proof: Looking at the explicit expression of Wigner's D-matrices ([3]) 

fc=max(0,2g) ^ ' ^ ^ 

we conclude that it is enough to have: |xn| < 1, \x 2 2\ < 1 and |a; 12^21 1 < 1, for the 
convergence of (jlj) for t = 1, because their exponents run up to infinity independent of 
each other. Moreover, if we require convergence in the expansions (flOj) and (fl2~l) for t = 1, 
then I det(X)| < 1 is needed too. ■ 

We shall see that Z and Z fulfil (19"B1) . but before we shall prove that 
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Proposition 6.2. For any matrix Z G ©4 we have that the squared norm of their rows 
is lesser than 1, that is: 

Nil! 2 + M 2 < 1, \z 2 l\ 2 + \Z 2 2? < I- (100) 



Proof: The positivity condition (|54|) says that 

det(.T - ZZ ] ) > <£> ku^i + ^i2^ 2 r < (1 - kn| 2 - ki 2 | 2 )(l - \z2i\ 2 - \z 2 2\ 2 ). (101) 

Hence, the last two factors must be either positive or negative. Supposing that both 
factors were negative would contradict ti(ZZ') < 2 in ( 1571) . Therefore, we conclude that 
both factors are positive. ■ 

Let us remind that, since Z and Z belong to D 4 , they must satisfy | det(Z)| < 1 and 
I det(Z)| < 1, as we saw in (l55l) and (|58l) . Now we are in condition to prove that: 

Proposition 6.3. The matrix X = Z^Z verifies the convergence conditions $98\) for every 
Z, Z G 4 and, therefore, the expansion §7ffl is well defined for Z = BD^ 1 . 
Proof: The conditions ( 11001) imply in particular that \zn\ < 1, \zizl < 1> \ z 2i\ < 1? \ z 22\ < 
1. Using this fact, the triangle inequality and taking into account that Z, Z G D verify 
( 11001) and the determinant restriction ( 1581) . we arrive to: 



\xn\ 


= l^llZH + ^12^21 


< 


. Z\ 1 Z\ 1 


+ 


| #12221 1 


< 


\z\\ 


+ 


1^12 1 


< 


\zn 


2 


+ 


\Z12 


2 


< 1 


\X 2 2\ 


= 521^12 + ^22^22 


< 


^21^12 


+ 


^22^22 


< 


\Z21 


+ 


^22 


< 


\z>2\ 


2 


+ 


\Z22 


2 


< 1 


\ X 12\ 


= \Zl\Z\2 + ^12^22 


< 


\ZwZ\2\ 


+ 


^12^22 


< 


\z n 


+ 


\zn\ 


< 


\z~11 


2 


+ 


\Z\2 


2 


< 1 


\x 2 l\ 


= \Z2\Z\\ + ^22^21 


< 


\z>2\Z\\\ 


+ 


-222^21 


< 


\Z21 


+ 


\Z22\ 


< 


\Z21 


2 


+ 


\z>22 


2 


< 1 



det(X) 



det(&Z)\ = I det(Z t )|| det(Z)| = det(Z t Z) 1/2 det(ZZ t ) 1/2 < 1, (102) 



which proves the convergence conditions ( |98i) .l 



7 Conclusions and Outlook 

We have constructed the CWT on the Cartan domain B 4 = U (2, 2)/U(2) 2 of the conformal 
group 50(4, 2) = SU(2, 2)/Z 4 in 3+1 dimensions. The manifold ©4 can be mapped one- 
to-one onto the future tube domain of the complex Minkowski space through a Cayley 
transformation, where we enjoy more physical intuition. This construction paves the way 
towards a new analysis tool of fields in complex Minkowski space-time with continuum 
mass spectrum in terms of conformal wavelets. It is traditional in Relativistic Particle 
Physics to analyze fields or signals (for instance, elementary particles) in Fourier (energy- 
momentum) space. However, like in music where there are no infinitely lasting sounds, 
particles are created and destroyed in nuclear reactions. A wavelet transform based on 
the conformal group provides a way to analyze wave packets localized in both: space and 
time. Important developments in this direction have also been done in [Ml [T5l fl~6] for 
electromagnetic (massless) signals. 
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In the way, we have stated and proved a A-extension ( TTUj) of the Schwinger's formula 
01]). This extension turns out to be a useful mathematical tool for us, specially as a 
generating function for the unitary-representation functions of S77(2,2), the derivation 
of the reproducing (Bergman) kernel of L^(04, g?z/a) and the proof of admissibility and 
tight frame conditions. The generalization of this theorem to matrices X of size N > 2 
follows similar guidelines and the particular details are discussed in the Appendix \K\ using 
the general SU(N) solid harmonics V^ a (X) of Louck [25]. This result could be of help 
in studying the discrete series (infinite-dimensional) representations of the non-compact 
pseudo-unitary groups SU(N,N). 

The next step should be the discretization problem. References [Ml EH] [36] give us the 
general guidelines to construct discrete (wavelet) frames on the sphere and the hyperboloid 
and [37] on the Poincare group. The conformal group is much more involved, though in 
principle the same scheme applies. 

Looking for further potential applications of the conformal wavelets constructed in 
this article, we think that they could be of use in analyzing renormalizability problems in 
relativistic quantum field theory. When describing space and time as a continuum, certain 
statistical and quantum mechanical constructions are ill defined. In order to define them 
properly, the continuum limit has to be taken carefully starting from a discrete approach. 
There is a collection of techniques used to take a continuum limit, usually referred as 
"renormalization rules" , which determine the relationship between parameters in the the- 
ory at large and small scales. Renormalization rules fail to define a finite quantum theory 
of Einstein's General Relativity, one of the main breakthroughs in Theoretical Physics. 
The replacement of classical (commutative) space-time by a quantum (non-commutative) 
space-time promises to restore finiteness to quantum gravity at high energies and small 
(Planck) scales, where geometry becomes also quantum (non-commutative) [35] • Confor- 
mal wavelets could also be here of fundamental importance as an analysis tool. 
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A Extended MacMahon-Schwinger's Master Theo- 
rem for Matrices of Size N > 2 

We have shown the utility of the Theorem 12.31 in dealing with unitary representations of 
SU(2, 2), in particular, in proving the admissibility condition 15.41 We would like to have a 
generalization of Theorem l2.3l for matrices of arbitrary size N, since it would be a valuable 
tool as a generating function for the unitary-representation functions of SU(N, N). 

The first step is to generalize Wigner D-matrices. This generalization has been done 
in the literature (see [25] and references therein) by the so-called SU(N) solid harmonics 
V^JX) defined as: 

T^{x) = ^m Yl ^ ( 103 ) 

A&M p N (a,/B) 

where the following space saving notations are employed: A is a N x N matrix in the non- 
negative integers a^] A\ = nfj=i a «r j X A = Yli j=% x ij J > a = ( a i> a 2, «at) is a sequence 
of N nonnegative integers that sum to p (i.e., a composition of N into p nonnegative 
parts), shortly a \- p; a\ = Ylf=i a t-'i M^-(a, p) denotes the set of all matrices A such that 
the entries in row i sum to on and those in column j sum to Pj, with a hp and P h p. 
Hence, T> p a JX) are homogeneous polynomials of degree p in the indeterminates Xij. 

The particular identification with Wigner's D-matrices for X of size iV = 2 is given 
by V{ uq2 {X) = V%(X) with p = 2j, a,P h 2j, a = (j + q u j - Ql ), P = (j + q 2 ,j - q 2 ). 
Matrices A G M%(a, P) can be then indexed by an integer k 

A^=( . k , i + <*- k ) (104) 
V J + <?2 - k k- q 1 -q 2 J 

with max(0, qx + q 2 ) < k < min(j + q%,j + q 2 ). 

The multiplication property OH]) and the transpositional symmetry (j7|) for Wigner 
matrices are still valid for SU (AT)-solid harmonics as: 

Y, VP ^ X ) V U Y ) = V U XY ) ( 105 ) 

a\-p 

and 

Kp(Y) = V%a(Y T ) (106) 

(see [25J for a combinatorial proof). 

Moreover, for general N x N matrices X, the determinant det(J— X) can be expanded 
in terms of sums of all principal q-th minors of X as 

N 

det(I -X) = d * det W' 

q=0 a\\-q 
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where the N- dimensional multi-index a = («i, (X2, ajv) is a partition of q with «j G 
{0, 1}, a fact that we now symbolize as a lh g; x = (in, x 22 , . . . , xjviv) and <9" = n»=i 
Let us define the sum of all principal (X — g)-th minors of X by 

T q (X) = ^det(X). 

att-N-q 

They are homogeneous polynomials of degree q = 0,1 ... ,N in the indeterminates x^. 
For example: T (X) = l,T 1 {X) = tr(X), . . . , T N {X) = det(X). Thus, det(7 - X) can be 
written in terms of these homogeneous polynomials as: 

N 

det(/-X) = ^(-l)%(X). (107) 

q=0 



Other possibility could be to use Waring's formulas 

To arrive at the A-extended MacMahon-Schwinger's Master Theorem (MSMT) for 
N x N matrices, we shall now proceed step by step from A = 2 to general A. Before, let 
us explicitly write down the generalization of the Theorem 12.11 to matrices of general size 
X 

Theorem A.l. (MSMT) The identity 

oo 

E tP E p -( x ) = det ( J - tx ) _1 ( 108 ) 

p=0 ahp 

holds for any N x X matrix X. 

The action of the operator Di on both sides of the Basic MacMahon-Schwinger's 
formula (11081) now gives: 

where we have defined T q (X) = (— l) q (q — l)T g (X). We can bring the numerator of the 
right-hand side of (I109p back to the left-hand side by using the expansion: 



1 _ I2q=2 T <l( tX ) p=0 q=2 

where we have used the following shorthand for 



oo N oo /sr^N \ 

n ^77^ = E(E f ^r = E J ; 2 13 f m 7 > 

p=0 q=2 7=0 ^ ^ ^ 



(110) 



T(xy = t 2 (x) 72 t 3 (x) 73 . ..T N (xy N , ( ^= 27j ) = ^3=^3 )• 



7 / 72! •••7at ! 
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Note that T{X) 1 are homogeneous polynomials of degree J2f =2 jlj in x ij- Inserting the 
expansion f II 1 j) in (11091) we conclude that 

oo oo ( /sr^N \ ^ 

E(f +DE H 27, j WEpl.(»o= di2) 

p=0 7=0 L v ' 7 J ahp v y 

This is the generalization of (1131) for general N. Let us proceed by applying D 2 on both 
sides of the identity (11121) : 

Y> + d E ( E ': 27) ) &> + 2 + E^«rEw = 2 2r^? • 

p=0 7=0 V ' 7 j=2 ahp V ; 

Using again (111 Oft , we have 

oo oo oo /s^N \ /ST^N i\ N 

E ^ E E F 1 : 2 7j J ( E ': 2 7< ) (p+ 2 + E i7,)f («)^ £ 

p=0 7=0 7 '=0 ^ I / \ I / j = 2 a hp 

' (113) 



det(J - tX) 

Rearranging series as in (|T5|) and making the change of (A/ - — l)-dimensional multi-index: 
a = 7 + 7', we obtain 

E ^ E IE p: 2 7J ) F^E 7J °) 0. + 2 + £ } 

p=0 (7=0 1 7=0 \ ' / \ I / j =2 ) 

*f( t xr^ Vaa (tx)= \ (in, 

a\-p ' 

Applying now D 3 on both sides of (11131) results: 

00 n 00 00 /sr^N \ /sr^N t\ N N 

£ ^ £ £ ( E ': 2 % ) F% 2 7 ') (p+2+E«)(p+ 3 +E «■* + ^» 

p=0 7=0 7'=0 ^ j=2 j=2 

1 — T (tX) 

xf (tX)^' J>L(«0 = 3 de ^ t t 4 

and using again f 1 1 1 1) we get: 

00 00 00 /\—\N \ /\~~^N /\ /t~\N //\ iV 



00 1 00 00 00 /v^jv \ /v^jv /\ /v^Jv //\ jv 

£ ^ £ £ £ ( 7 2 7 ( y 7i E ' y ? % >+ 2 + E«) 

p=0 ' 7=0 7'=0 7"=0 V I / \ I / \ 1 / j =2 
N ~ 1 

x(p + 3 + £j(7i + T ;))T(tX)^" J>L(**) = det(/ _ tx)4 - ( 115 ) 
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Rearranging series and making the change a = 7 +7' +7", we can recast the last expression 
as: 

3 Z^]2^^l 7 A 7' / V o- - 7 - 7' 

p=0 (7=0 L 7=0 7'=0 V ' 7 V ' 7 V ' ' 

- ^ 1 ~ 1 

x (p + 2 + £ J7i)(p + 3 + E^ + 7;)) T(t*)" E^tX) = detf/ _ tx u - 

j=2 j=2 J ahp 1 ^ 



If we repeat the process (A — 4) more times, then we arrive at the following identity: 

00 , 00 00 00 /x—^N (A— 2)\ A— 3 /k—kN (k) 



p=0 v ' 7=0 7 '=0 7 (A-2) =0 v ' ' fe=0 

AT k 

:(p + k + 2 + J E7f )f (tX)^'+-^ (A - 2) ^PL(tX) 
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, n dot (7 - / .V)- 

j =2 1=0 ahp v 7 



Making once more the change a — 7 + 7' + • — h 7^ 2 -*, we can write: 

E B E E - a^F- 

p=0 (7=0 ahp v y 

where we have defined the following coefficients: 

a (7-7 <T- 7 7 (A_4) /ST N (rr _ V A " 3 -v W ' 

_ 1 \^ V- \ ^ ll^j=2\ U i Z^i=0 'i . 

P ' CT ~ (A-2)!^^' ^ V a-V A - 3 7« 

v ; 7=0 7 '=0 7 (A-3) =0 v u Z^,i=0 I 

A-3 /s~^N (fc)\ N k 

fc=0 V ' 7 J=2 1=0 

In order to account for the particular coefficients C 2 i(T and Cp a , given inside curly brackets 
in (I112p and HI 14)) . we must understand in (11161) that: 1) summations on 7® with k < 
are absent, 2) empty or miliary sums are zero and 3) empty or miliary products are 1, as 
customary. 

Summarizing, we can enunciate the following: 

Theorem A. 2. (A-extended MSMT) For every A £ N, A > 2 and every N x N matrix 
X , the following identity holds: 



E xett E^-^c-^f (x)- y, v l {x) = det(/ - tx) " A ' (117) 



p + 1 

A- 1 

p=0 (7=0 ahp 

with C x pa given by HIE) and f(X) a by 17771) . 
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The expression (11171) generalizes ( TTUj) for matrices X of arbitrary size N. In fact, for 
N = 2, the coefficient (11161) reduces to: 

MV-7 2 )( A - A + -r 2 > ™ 

which agrees with ( FlO]) . We have also been able to find simplifications of a in the follow- 
ing cases (we take the binomials in the generalized sense (£j = n(n — 1) . . . (n — m + l)/m\ 
to account for fractional n): 

i) For 2 < A < 5, the coefficients (11161) are given by: 



'p,a 



( a - 2 + p + 1 ^;f =2 fc(T^ i 



A-4 



' ' fc=3 J 



ii) For iV = 3 and A > 2, the coefficients (I116p can be obtained trough the expression: 



x 



&2 + cr 3 \ /A - 2 + (T 2 + 0"3 
a )\ A-2 

i=l ^ ' j=l ^ 



where we have defined £ = Odd(A), that is, £ = when A is even and £ = 1 when 
odd. See that ffTT9l reduces to ( FTTH]) for <x 3 = 0. 



B Continuous Wavelet Transform on a Manifold: a 
Brief 

The usual CWT on the real line M. is derived from the natural unitary representation of 
the affine group G = SIM{1) in the space of finite energy signals L 2 {R, dx) (see Section 0] 
for a reminder). The same scheme applies to the CWT on a general manifold X, subject 
to the transitive action, x — > gx, g G G, x G X, of some group of transformations G 
which contains dilations and motions on X. If the measure dv(x) in X is G-invariant (i.e. 
du(gx) = du(x)), then the natural left action of G on L 2 (\, dv) given by: 

[U(g) ( p}(x) = ( p(g- 1 x), geG,^eL 2 (X,du), (120) 

defines a unitary representation, that is: 

(U{g)y\U{g)ct>) = {y\cj>) = [ ^x)<j>(x)du(x). 

Jx 
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When dv is not strictly invariant (i.e. dv(gx) = j\4(g,x)dv(x)), we have to introduce a 
multiplier ( Radon- Nikodym derivative) 

[U(g)4>](x) = M^xf'^g-^x), geG,4>e L 2 (X,dv), (121) 

in order to keep unitarity. The fact that U(g)U(g') = U(gg') (i.e. U is a representation of 
G) implies cohomology conditions for multipliers, that is: 

M(gg',x) =M(g,x)M(g',g- 1 x). (122) 

Consider now the space L 2 (G,dfi) of square-integrable complex functions ^ on G, where 
dfi(g) = dfi(g'g), Wg' G G, stands for the left-invariant Haar measure, which defines the 
scalar product 

(*|$) = / ^g)$(g)dfi(g). 

JG 

A non-zero function ip G L 2 (X, dv) is called admissible (or a fiducial vector) if $?(g) = 
(U(g)t/)\i[>) G L 2 (G,d/i), that is, if 

c^= f ^g)^{g)df,(g) = [ \{U{g)^)\ 2 d^g) < oo. (123) 

JG JG 

A unitary representation for which admissible vector exists is called square integrable. 
For a square integrable representation, besides Eq. (1123!) the following property holds 
(seeP): 

\(U(g)m\ 2 M9) < oo,V0 G L 2 (X,du) . (124) 

Let us assume that the representation U is irreducible, and that there exists a function 
ip admissible, then a system of coherent states (CS) of I/ 2 (X, dv) associated to (or indexed 
by) G is defined as the set of functions in the orbit of ip under G 

^ g =U(g)^P, geG. 

There are representations without admissible vectors, since the integration with re- 
spect to some subgroup diverges. In this case, or even for convenience when admissible 
vectors exist, we can restrict ourselves to a suitable homogeneous space Q = G/H, for 
some closed subgroup H . Then, the non-zero function %p is said to be admissible mod(i/, a) 
(with er : Q — > G a given section) and the representation U square integrable mod(if, a), 
if the condition 

/ \{U(a{q))i)\<P)\ 2 dfi(q) < oo, V0 G L 2 (X,du) (125) 
JQ 

holds, where dp, is a measure on Q "projected" from the left- invariant measure d\i on the 
whole G (see [39] for more details on this projection procedure). Note that this more 
general definition of square integrability includes the previous one for the trivial subgroup 
H = {e} and a the identity function. The notions of square integrability and admissibility 
mod(-ff, a) were introduced in [40] (see also [4]). 
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The coherent states indexed by Q are defined as ip a (q) = U{a{q))ip,q G Q, and they 
form an overcomplete set in L 2 (X, du). 

The condition ( 11251) could also be written as an "expectation value" 

0< / \{U(a(q))^\<P)\ 2 djl(q) = {(j)\A<r\<f) < °o, V0 G L 2 (X, di/), (126) 

where ,4 CT = fQ\4>a(q))(4>a(q)\d(i(q) is a positive, bounded, invertible operator. If the 
operator A' 1 is also bounded, then the set F a = {\if>a( q )), q G Q} is called a frame, and a 
frame if At is a positive multiple of the identity, A a = cl,c> 0. 
To avoid domain problems in the following, let us assume that ip generates a frame 
(i.e., that A" 1 is bounded). The Coherent State map is defined as the linear map 

%: L 2 (X,dv) — ► L 2 (Q,djl) 

<j> i — > %4> = 

with $w,(g) = it s range L\{Q,dfi) = %b(L 2 CK, du)) is complete with respect to 



the scalar product ($|$')^, = (^T^pA^T^ 1 $') n and 7^ is unitary from L 2 (X,du) onto 
L 2 (Q,<i/2). Thus, the inverse map 7^ _1 yields the reconstruction formula 

4> = X p - 1 %= I Q*(q)AzH*{ q )dp.(q), $4,eLl(Q,dp), (128) 

J Q 

which expands in terms of coherent states (wavelets) A~ 1 ipcr(q) with wavelet coefficients 
^^{q) = [T^^il) ■ These formulas acquire a simpler form when A a is a multiple of the 
identity, as it is precisely the case considered in this article. 

C Ort honor mality of Homogeneous Polynomials 

In order to prove the orthonormality relations ( |66l) . we shall adopt the following decom- 
position for a matrix Z G ©4 

Z = UxSUl 

where U lf2 G U{2)/U{1) 2 [as in with ft = ft = 0] and 3 = diag(£ x , 6},£i, 2 G D x . 
This parametrization ensures that Z G D 4 ; in fact 

I — ZZ^ = Ui(I — EE^)U{ > (129) 

since the eigenvalues are 1 — | £1,2 1 2 > 0. 

Let us perform this change of variables in the invariant measure (1M|) of ^(D^di/x). 
On the one hand, the Lebesgue measure on C 4 can be written as: 

\dZ\ = J|d£i||d&l<te(tfi)<fe(t/-2), 
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with ds(Uift) defined in (162]) and J = |(|£i| 2 — K2I 2 ) 2 is the Jacobian determinant. The 
Lebesgue measures {d^i^l and lebi^l will be written in polar coordinates = Pk^ ldk and 
Zk = rke iak , k = 1, 2. On the other hand, the weight factor in (|64l) adopts the form 

det(J - ZZ^ = ((1 - pj)(l - pl)) x - 4 = fi(p), 

so that the invariant measure of L^(D^,dvx) reads: 

du x (Z,tf) = cxJipMp^Wd&ldsiU^dsiUi) (130) 
= cxJ(p)Q(p)pidpid6 1 p 2 dp2d8 2 (l + r\)~ 2 ridridai(l + rl)~ 2 r 2 dr 2 da 2 , 

with < p 12 < 1, < r x 2 < 00, < 0i 2 < 27r, < 2 < 27T. Let us call 



2 j + 1 /m + A - 2\ /m + 2j + A - 1 



Mj ' m ' » A-l V A - 2 y V A - 2 
the normalization constants of the basis functions fl65l) . We want to evaluate: 

feKS> = / ^)det(Z)»»x4^(Z) det(Z)-'Pj; 9 , (Z). (131) 

Using determinant properties, Wigner's X>- matrix properties (j6]) and (j7]), and the fact 
that det(L r i !2 ) = 1 and S is diagonal, the previous expression can be restated as: 

<W' m tn j '' m '\ $ 3 



ni 



9= -j <?'=-y 1 



EE C J Jn|d6 L ||de 2 |2? J M (S)I>j; f ,,(S)det(Srdet(s; 



X 



da(CTi)^(l/i)2^ >9 ,(l7i) / ds{U 2 )V^ q {U 2 )% . (U 2 ) (132) 



Let us start evaluating the first integral. For the diagonal matrix S we have that 
W iig2 (^) = 5 gi , g2 ^' +9l er <?1 ,sothat 



— p> + 3"+<l+<l'+m+m' pj+j'-q-q'+m+m' e i(j'-j+q'-q+m'-m)6i e i(j'-j+q-q'+m'-m)d2 (133) 



Integrating out angular variables gives the restrictions 

r r V3 qq (s)v^ ( =■) det(^r det(sr' de.de, 

Jo Jo 

— 47r Oq,q'Oj+m,j>+m>Pl P 2 

Integrating the radial part: 

4tt 2 c a t C J(p)n(p)pf j+q+m) p?- q+m) p 1 dp lP2 dp 2 



jo 



(j + m) 2 + (j + m + 2g 2 + 1)A - 5g 2 - 1 = , 

7T 2 (A - l){3 +m+ 1 +X - 1 ^ (3+m-q+\-l^ ~ '^j+m 
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and putting all together in f 1 1 3 2 j) we have: 



(, j,m ,J'> m '\ min{j,i'} 

" j+m,j' '+m' / j '^j 



J V j, rrw » j'.m' . 

?=-min{3j'} 

x / ds(l/ 1 )x4 >fl TO2^ | ff (l7i) / ds{U 2 )TP (U 2 )VZ (pi) (134) 

The last two integrals are easily computable. Actually they are a particular case of the 
orthogonality properties of Wigner's D-matrices. More explicitly: 

ds(U)W a AU)V\ (U)= / , " X>2 a {UW, (U)=5 ii ,6 <na .——. 
Going back to f)134p it results: 

Finally, taking into account the combinatorial identity: 

2j + 1 



J2 (A - 1)^ 



9=-J 



/m+A-2\ /m+2j+A-l\ 
I A-2 / I A-2 / 



and the explicit expression of the normalization constants A/}, m , we arrive at the orthonor- 
mality relations ( )66l) . 
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